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We consider a gedanken experiment with a beam of atoms in their ground state that are accelerated 
through a single-mode microwave cavity. We show that taking into account of the "counter-rotating" 
terms in the interaction Hamiltonian leads to the excitation of an atom with simultaneous emission 
of a photon into a field mode. In the case of a slow switching on of the interaction, the ratio of 
emission and absorption probabilities is exponentially small and is described by the Unruh factor. 
In the opposite case of sharp cavity boundaries the above ratio is much greater and radiation is 
produced with an intensity which can exceed the intensity of Unruh acceleration radiation in free 
space by many orders of magnitude. In both cases real photons are produced, contrary to the opinion 
that a uniformly accelerated atom does not radiate. The cavity field at steady state is described by a 
thermal density matrix. However, under some conditions laser gain is possible. We present a detailed 
discussion of how the acceleration of atoms affects the generated cavity field in difi'erent situations, 
progressing from a simple physical picture of Unruh radiation to more complicated situations. 

I. INTRODUCTION 

Intriguing properties of vacuum as viewed by accelerated observers have been the subject of intense investigation 
for almost three decades. One of the most remarkable results is the so-called Fulling-Unruh effect predicted and 
analyzed by Davies, Fulling, Unruh and DeWitt 1], and others In essence, it was shown that ground state atoms, 
accelerated through vacuum, are promoted to an excited state just as if they were in contact with a blackbody thermal 
field. These studies predict that a (two- level) ground state atom, having transition frequency w, and experiencing 
a constant acceleration a, will be excited to its upper level with a probability governed by the Boltzmann factor 
exp(— 27rci;/a), where a = a/c, c is the speed of light in vacuum. Unfortunately, even for very large acceleration 
"frequency" a « 10^ Hz 3], and microwave frequency oj « 10^" Hz 4], this factor is exponentially small, ~ 10^^""; 
and is not of experimental interest. 

To begin with, we note that the physical picture of the Unruh effect in any setting is quite straightforward. In 
particular, it is the counter-rotating terms in the interaction Hamiltonian that describe the process of an excitation 
of an atom with simultaneous emission of a photon (see Fig. 1). 

Furthermore, it is known that many effects related to the interaction of atoms with the electromagnetic field acquire 
new features or are enhanced in the cavity quantum electrodynamics (QED) setting, when the atoms are injected into 
a high-Q cavity. 

Thus we were motivated to study a simple gedanken experiment with a beam of atoms in their ground state that 
are accelerated through a single- mode microwave cavity 5]; see Fig. 1. This model is sufficiently simple so that we are 
able not only to find the probabilities of a photon emission and absorption by ground state atoms, but also to solve 
the density matrix equation for the photons in the cavity mode interacting with a beam of atoms and to analyze its 
steady state solutions. 

We found that the radiation is thermal (in the typical case) and the effective "Boltzmann factor" can be much 
larger than the above exponentially small value Q. In particular, for the above example it is given by a/2TTUj, which is 
of order ~ 10"'^. Hence, it is many orders of magnitude larger than that for the usual Unruh effect and is potentially 
observable. 

We find that the mechanism of simultaneous excitation of both field and atom in a cavity is the same as for the 
Unruh effect in free space. We show that in both cases it is the nonadiabatic transition due to the counter-rotating term 
a^cF^ in the interaction Hamiltonian The reason for an enhanced excitation in the cavity is the relatively large 
amplitude for a quantum transition |6, 0) — > |a, 1) between the dressed atomic states due to the sudden nonadiabatic 
switching on of the interaction at the cavity boundaries, whereas for the Unruh effect in free space the emission is 
exponentially small due to a slow switching on. In both cases nonadiabatic effects, however small they are, play a 
critical role: there is quite a real emission of a photon accompanied by the excitation of an atom - not just dressing of 
the ground state of an atom as a result of interaction. As was shown in |5jj| , when the cavity boundaries are removed, 
our expressions yield the usual Unruh result. 



2 



(a) 



b b 



b b 



(b) 
V 

V 



(c) 



AW 



0) = 



E-E, 




FIG. 1: (a) Resonant absorption or emission: an atom is excited (deexcited) as it simultaneously absorbs (emits) a photon, 
(b) Counter-resonant absorption or emission processes that are usually neglected in the "rotating wave approximation": an 
atom is excited (deexcited) as it simultaneously emits (absorbs) a photon, (c) The energy for counter-resonant emission is 
drawn from work done by a force accelerating an atom, (d) Atoms or ions in their ground state \b) are accelerated through a 
single-mode microwave or optical cavity. 



Other processes where counter-rotating terms play a crucial role include, e.g. parametric resonance and the anoma- 
lous Doppler effect. The similarity between the Unruh effect and the anomalous Doppler effect, in which an atom 
moving faster than the velocity of light in a medium is also excited to an upper state while simultaneously emitting a 
photon, has been emphasized in 6]. In both cases the energy for the excitation of an atom and emission of a photon 
is taken from the work done by the force supporting the motion of an atom along the given trajectory. However, for 
the Doppler effect there are no time-changing parameters: the excitation occurs due to the existence of the Cerenkov 
resonance. 

In section II we formulate the model and write down the Hamiltonian. In section III the transition probabilities for 
emission are shown to yield a simple physical picture of the Unruh radiation. Next, the master equation is derived 
and the steady-state solution for the photon density matrix is derived and analyzed. In sections IV- VI we analyze the 
mechanism of emission and absorption by the accelerated atoms and its relation to the standard Unruh effect, the 
emission and absorption probabilities are calculated. The resulting integrals are evaluated by the method of stationary 
phase in all physically interesting asymptotic limits. The possibility of amplification and laser action is discussed. 
The interpretation of the results and discussion are presented in section VII. 



II. THE MODEL 



We start from writing Hamiltonian for the system consisting of a two-level atom interacting with the electromagnetic 
field: 



H = Ha + Hf + V. 



(1) 
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Here Ha = hwa^ is the Hamiltonian for an atom with ground and excited states \b > and |a >, respectively, separated 
by energy difference Ea — Eb = hco, = l/2(|a >< a\ — \h >< b\) the Pauh matrix, and Hf = J2k ^^^fcO^Ofc is the 
field Hamiltonian, where and ak are photon creation and annihilation operators and the k-summation is taken 
over the electromagnetic modes of a cavity or a free space, depending on the formulation of the problem. 

The atom-field interaction Hamiltonian in the atomic rest frame can be written in the dipole approximation as 



V = hS^gk (flfc + flfe) Uk{z{T)) (a + (7+) 



k 



(2) 



Here Mfc(r) form a set of orthogonal, normalized functions, a'^ = \a >< b\ and a = \b >< a\ are the atomic raising 
and lowering operators, and gk = fJ.E'-y./h is the atom-field coupling frequency, which depends on the atomic dipole 
moment fi and the electric field amplitude E'y. in the rest frame of an atom, evaluated on the trajectory z(t) of an 
atom as a function of proper time r. 

In the interaction representation, the master equation for the density operator p can be written as 



dp 
in— — 
dr 



V{t)p 



(3) 



where the interaction Hamiltonian in this representation can be obtained by replacing Ofc exp(— z;/fei(r)), — > 

exp(wfet(T)), a — *■ (Texp(— iwr), and ct"*" — *■ exp(zcijr), where t is the time in the incrtial laboratory frame. Note 
that we do not intend to use rotating wave approximation since the excitation of an atom from the ground state with 
simultaneous emission of a photon is described by counter-rotating terms of the type d'^a^ . 

Our main goal will be to solve master equation Eq. ^ and analyze its solution in various limits. However, we can 
get some insight into the physics of an accelerated atom-field interaction by solving first a simpler problem, namely, 
calculating photon emission and absorption probabilities within the standard first-order perturbation theory. 

III. TRANSITION PROBABILITIES, MASTER EQUATION, AND THE PHOTON DENSITY MATRIX 

A. Probabilities of photon emission and absorption 

Consider an atom entering the cavity at a proper time r^, at which moment the interaction with a cavity mode 
is assumed to be turned on. If the interaction is weak enough, the state vector of the system atom-|-field at any 
subsequent time t can be found using the first-order perturbation theory: 



The probability of transition |'0(Ti)) is therefore given by 

1 



P 



ft2 



{^{T)\V{T')\tP{n))dT' 



(4) 



(5) 



In particular, if an atom was initially in its ground state \b >, the probability of excitation of an atom with 
simultaneous photon emission into the kth mode can be calculated as 



1 



gk{r'){lk\a+\Ok){a\a+\b) dr' 



(6) 



The probability of photon absorption from the fcth mode by a ground-state atom, when there is only photon in this 
mode, is given by 



P(Ofc,a) = l 



gk{T'){Ok\ak\lk){a\a+\b) dr' 



Now consider a uniformly accelerated atom moving along the trajectory |l2l | 

1 c 

t(r) — to -\ sinh(aT), z(t) — — [cosh(aT) — 1] . 

a a 



(7) 



(8) 
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where = t{T = 0) is the moment of time in the laboratory (cavity) frame when the atom starts its acceleration. 
Expanding the electromagnetic field in terms of running waves with wave vectors k, fc^ = k • v/w, the atom- field 
interaction Hamiltonian in the atomic frame is given by 

V{t) = J2 ^9kM [afee-^'^*(^)+*'=-^(^) + h.c.j [ae-''^^ + h.c] . (9) 

For the dipole moment fi oriented in x-direction, the corresponding electric field amplitude in the atomic frame E'^ 
is related to the x-component of the electric field amplitude in the lab frame as E'^ = a/(c — u)/(c + v)Ex. Since 
V — c tanh(Q!r) for a uniformly accelerated particle, we have E'^ = e~°''^ E^ and gki^) = f/fee"""^. 

For simplicity, consider the case of a single-mode cavity and co-propagating atom and field: kz — |k| — v/c. 
Substituting Eqs. ijHJl, (HTJ into Eqs. ©, Q, we obtain 

P(0fc,a)=g2|j^(^)|2. p{i^^a)=g^\h{u;)\\ (10) 
where the absorption and emission amplitudes are given, respectively, by 

/a(w) = exp [i^ (e""^ - 1) + iwr - ar] dr, , , 

/,(c.)=/;(-^), ^''^ 

and it is assumed that an atom enters the cavity at r = and exits cavity at r = Tg. Hereafter we skip the index k 
in the coupling constant gk for the single- mode cavity. For a counter-propagating wave one needs to replace a —^ —a 
in Eq. l(TT|l. 

We also note that the amplitude of the process of photon emission by an excited atom is = I*. Therefore, the 
probability of this process is P(l/c, h) = P^Ok: a)- 

If we remove the cavity walls to infinity by letting — > — cx) and Te — > cx), the integrals in Eq. Hll|l are reduced to 
gamma-functions by making the substitution x = ^e"""^: 

/..e = ie-'^^ (^) r e-x^^^ d. = le-'^^ (^) %±i^r f 1 T ^) , (12) 
ly yiy / Jq v \v / \ a / 

where the upper and lower signs correspond to the absorption and emission amplitudes respectively. Using the equality 



we arrive at 



oe~ — r (14) 



27TUJ 



sinh(^) 

and 

Note the familiar Planck factor (^e'^B^^ — 1^ in the emission probability (|15l) with temperature equal to the Unruh 

temperature T„ = 2^kBC - ratio of emission to absorption probabilities is also in agreement with Unruh result in 
free space: 



P(Ofe,a) 



g-^TTc^/a^ (16) 



It is interesting that the simple one-dimensional model of ref. and the present discussion contain the basic physics 
of the Unruh effect even in the free-space limit. Strictly speaking, in that limit, the interaction of the atom with all 
wavevectors k of the vacuum field needs to be taken into account. In Sec. IV and V we provide a detailed interpretation 
of this result and also generalize the treatment to include the modes with different orientations of k. 
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B. Master equation for the density matrix 



Consider again uniformly accelerated atoms moving along a trajectory Q and interacting with a single-mode field 
in a cavity. Our goal is to find the solution to the master equation Eq. Q within the perturbation theory following 
the approach described in [Tol lllj | . In particular, we will find the steady state number of photons in a cavity mode 
as a result of interaction with a beam of atoms. 

Atom in free space. 

First, we derive the results for a single atom in free space within the approach of the quantum theory of the laser 
[Ullll. Equation © can be rewritten in the form convenient to apply perturbation theory expansion: 



dp 
d^ 



V,p 



V{t), F(r'),p(r') 



dr'. 



(17) 



In Markov approximation, assuming weak interaction with many field modes we decompose the density matrix as 
the product of its atomic and field parts as p(t) ~ Patomij) ® pfid{Q). Tracing over the field degrees of freedom we 
obtain from Eq. (^Tjl 



dpa 



dr 



1 



Tr. 



fld 



V{t)V{t')p 

atom Patoin 



-V{T)patom{T)pfld{m{r') " V {t') Patomir) p fld{m i^) 



dr'. 



(18) 



Let us first consider an atom at rest, when t — t. In the interaction representation, using Eq. and the replacement 
described after Eq. we arrive at the following equation for population of state \a>: 



dpa 



dt 



2 pt 



dt' 



Eld,{t)dl{t'))e^-'^'-'"^ +{t^ t'} ] Pa, 
{J2Eld+it')d,it))e^-^'-''^+{t^t'} 



Pbb 



(19) 



where the field operators in angular brackets can be expressed via average number of photons in the fcth mode as 



(a+(t')afc(t)) =nfce— 
{ak{t)d+{t')) = (l + nfc)e— 

After performing integration, we finally get Eq. H19|) in the form 



dpaa 

dt 



= -const X ^ [(1 + nk)paa{t) + nkPbbit)] 



(20) 



(21) 



Its steady-state solution when fife is a thermal field is Paa/ Pbb = . 

For an accelerated atom in Minkowski vacuum one can obtain a familiar Unruh result. One can also obtain a more 
general result for an atom accelerated through a thermal (not vacuum) background electromagnetic field with photon 
distribution rife: 



dpaa 

dt 



^t' ^El{nu\{dk + d+)r{dk + d+)'^\nk) 

k 

/ dr' |nfce*''''(*'*')~''''("^(*)~"^(*')) -|- (1 + nfe)e^*''''(*~*')+*''(''(*)""^(*'))| 



c.c. 



(22) 



where t = t{T),t' = ^(t'). 

Next, we proceed following the method described in e.g. P. Milonni or J. Audretsch and R. Miiller 2]. Namely, 
we assume that the frequency Vk has a small imaginary part, substitute the equations for a uniformly accelerated 
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trajectory t(T), z(r) and perform a summatfon over k which leads to 



dT 



dr' 



nkn 



1 + 



(sinh(a(r — t')/c + iae/ c)Y (sinh(a(T — t')/c — iae/c)Y 



(23) 



where fco = lo/c. Then we represent functions l/(sinhx)^ as infinite series X]p {a{T~T') /c~-Kip±ieY ' ^-'^P^^^ time 
integration over infinite hmits, evaluate the integrals by method of residues and finally arrive at 



dpaa 
dT 



= -/3 [nxnA + [riT + l)inA + 1)] Paair) + /3 [nrinA + 1) + [riT + Pbbir), 



where 



e = — 1 



-; nr 



QkT — 1 



(24) 



(25) 



and /3 is a constant which is unimportant for a steady-state distribution of populations. 

Equation (|24|l allows one to find steady-state atomic populations for a general case of an atom accelerated through 
a thermal field background. 

Beam of atoms accelerated through a single-mode cavity. 

In this part our goal is to evaluate the steady-state number of photons in a cavity mode as a result of interaction 
with a beam of atoms. 

As in the quantum theory of the laser [Tol [T]| , the (microscopic) change in the density matrix of a cavity mode due 
to any one atom, Sp^, is small. The (macroscopic) change due to AN atoms is then Ap — Sp^ — ANSp. Writing 
AA^ = rAt, where r is the atomic injection rate, we have a coarse grained equation of motion: Ap/ At — r5p. The 
change 5p^ due to an atom injected at time in the atomic rest frame is 



- -w 



Tr, 



(26) 



V{t'), T/(r"),p'^'°"(r.)®p(i(r,)) 



dT'dr", 



where tratom denotes the trace over atom states. The time r is the atomic proper time, i.e., the time measured by an 
observer riding along with the atom. For simplicity, consider again the case of the co-propagating atom and field and 
the interaction Hamiltonian given by © . 

In the case of random injection times, the equation of motion for the density matrix of the field is 



dpn,n/dt = -i?2 [{n + l)pn,n - npn-l^n-l] 
[npnji - (n + l)pri+l,,i+l] , 

where R12 are defined in the following. If i?i > i?2, there is a steady state solution which is thermal ^ 



(27) 



-hun/kBT , 



.^—hu/ki 



(28a) 



n 

where an effective temperature of the field in the cavity is T = hv /kB^T^[Ri/ R2]- Thus, spontaneous emission of 
randomly injected ground state atoms in the cavity results in thermal statistics of the mode excitation. Note, that the 
thermal statistics of the atomic excitation in the standard Unruh effect in free space is due to spontaneous emission 
into a vacuum field reservoir with a continuous spectrum of modes. 

Absorption and emission coefficients i?i_2 = ^"1.9-^1, 2 P are determined by the amplitudes = ~j[ It'^'^ Vi,2dT of 

the matrix elements Vi = (a, 0|F|6, 1) and V2 = (a, llV^I^, 0) of the interaction Hamiltonian (jSJ, respectively, and their 
explicit form is given by Eq. Hll|) . Using the results of the previous section, we get the same result that in the limit 
v,uj^ a the emission/absorption ratio is R2/R1 ~ a/(27rti;), which is an enhancement by many orders of magnitude 
as compared to the exponentially small value R2/R1 = exp(— 27raj/a). 
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IV. EMISSION AND ABSORPTION OF RADIATION BY GROUND-STATE ATOMS 



After the substitution of variables x — ^e "'^ , the absorption and emission ampHtudes can be expressed via 
incomplete gamma- functions: 



e±5S 



-a(Te 



-•))-r(^,«) 



(29) 



where ^ = 1 =F i^, u = — J^e"""^' , and r(^, u) = e~^x^~^dx is the incomplete gamma function. 

In principle, expressions H29|l can be fully analyzed because the properties and asymptotic behavior of incomplete 
gamma-functions are well known. Some representative graphs of the emission and absorption amplitudes as functions 
of the field frequency will be shown below in Figs. 2,3,4. However, it is more instructive and transparent to directly 
calculate the asymptotic of the integral by applying integration by parts and the method of stationary phase. 

In particular, we consider the most realistic case iy,Lu ^ a and apply the stationary phase method that can be 
summarized as 



.5 

/ F{T)e''^f^^'>dT = B + S, 

J a 



where 



B 



F{r)e 



lAf'ir) 



N 

E 



-1 d 



\ f'{T) dT 



F{T)e 



iAf(T) 



fir) 



o{A 



-N\ 



(30) 



(31) 



is the contribution from integration boundaries obtained by integration by parts. 



I 27ri 



(32) 



is the contribution from a stationary point Ts such that /'{ts) = 0, f"{Ts) ^ 0, obtained by expanding fir) in Taylor 
series around Tg. It is assumed that A ^ \. We will also consider separately the case when the stationary point 
approaches one of the integration boundaries; see Eq. (|38|l below. 

Suppose for definiteness that v > uj and = 0. When v — lo > ^/acJ, the stationary point = ^ log ^ of the 
absorption integral la in Eq. is within the integration limits and far enough from the boundaries. Therefore, la 
can be evaluated as a sum of the boundary contribution 



/(&) 



exp 



L a V 



1) -|- iujTe 



1 



-lue 



i{v — Lv) 



and the contribution from the stationary point 



27r to 




auj\ V 



— exp 



V TT 

I- log - +«- 

a 4 



(33) 



(34) 



It is clearly seen that the contribution from the stationary point dominates in the absorption integral la- The same 
result can be of course obtained directly from Eq. (|29|l after moving the integration boundaries to ±cxi and considering 
the resulting expression 



/i(c.) 



V \ v J \ a J 



(35) 



in the asymptotic limit of a large complex argument of the gamma-function. 

The emission integral /g, which originates from the counter-rotating term cx a^&^ in the interaction Hamiltonian, 
does not have a stationary point within the integration limits. Therefore, its value is solely determined by the boundary 
contribution, and /e(i^) ^ Ia^\—uj). If we further assume long enough interaction time, ar^ 3> 1, the second term on 
the right-hand side of is much greater than the first term, and we obtain 



(36) 
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which is equal to 2^ for v ^ uj. 

Exactly at resonance, v — uj = ^auj, the stationary point coincides with the lower integration limit t — Q. In this 
case one can show that the main contribution again comes from the stationary point, and the value of the integral 
lie is two times smaller than (I34|l . The resulting ratio of probabilities is equal to 

P(Ofe,a) 2TIU- 

The above analysis can be readily generalized for an arbitrary value of ^^^=- In this case the stationary phase 
method gives an additional term in the integrals /i_2 that contains the error function crf[z] of detuning: 



/ /27rci; ( u — V uj v n 

- 1/1 r:^exp i + i-\og - + i- 

\ \auj\ ziy \ a a ui 4 



1 + erf ^ e 



(38) 



The function gives the spectral profile of the absorption line. 

The ratio H36(l or 1)3 7|l is surprisingly large; in fact, it is exponentially larger than the value 

one would expect to obtain on the basis of studies of the Unruh effect. Here 

ha 

ksTu = — (39) 

is Unruh temperature. In our case the effective temperature of radiation in the vacuum state of the cavity mode is 
determined from 

a 

e ^bt = — ^ 



which gives 



keT = ^ . (40) 

l0g222£ ^ > 



The reason for such a large effective temperature is apparently the sudden turn-on of the interaction of an atom 
with a cavity mode. If we eliminate the nonadiabatic switching effect by letting Ti —00 and Te — » 00, the integrals 
in Eq. (|^ are reduced to 



Using the equality 



we arrive at the Unruh-type result 



Ia,e = -e-'^ (-) e±^r It — 



V. ANGULAR DEPENDENCE OF EMISSION/ ABSORPTION PROBABILITIES 

The above conclusion does not depend on our assumption of interaction with a single co-propagating cavity mode 
and can be generalized for the case of an electromagnetic mode with an arbitrary k-vector. Similarly to Sec. Ill, 
we calculate the probability P(lk, a) of excitation of an atom with simultaneous photon emission into the kth mode 
assuming that the field was initially in the vacuum state. Then we calculate The probability P(Ok, a) of photon 
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absorption from the fcth mode by a ground-state atom, when there is only photon in this mode. The arguments of 
the P functions denote the final state of the field and atom. The ratio of these probabilities is given by 

P(Ok,a) /kH ' ^ ^ 

where 

Ik{Lu) ^ J — explivt^r) — ikzz{T) — iujT ~ ar] dr, (44) 

where A; = |k| = i//c. The probability of emission by an atom into all electromagnetic modes is proportional to J /k cPk. 
We will be interested in evaluating the ratio (|^ . Using equations for the trajectory of a uniformly accelerated atom, 
we arrive at 



/ exp i— I sinh QT ^coshar ) 

Jn I a \ k J 



lUJT — ar 



dr. (45) 



As in the case of co-propagating mode, the above integral can be calculated exactly in the infinite limits and can 
be evaluated approximately by the method of stationary phase in finite limits. 

For the infinite integration limits, Ti,Te oo, it was shown in that it is convenient to change the integration 
variable to 

l3 = ar ~ rj, 

where tanhry — k^/k. Then the integral in l|45() can be written as 

oo 

si„h/3-«/3-i^„-„ ^ 2e-''f "-"-^i^^(K^), (46) 

-OO 

where k± — k^c/a^ ^ = 1 + iou /a^ and K^{k±_) is McDonald function. Using the above result in (|43|1 . we obtain 

I 1 2 

f(lk,a) ^ |j^l-ic^/a(K±)| , . 

no.,a) ' " |K,+../.(.,)r ^ ^ 

which is "almost" Unruh factor in the limit Lo/a ^ 1, since K^p{x) — Kp{x). The extra factor of 1 in f is due to 
the fact that we are dealing with photons that have spin 1. This introduced an additional term in the integral as a 
result of Lorentz transformation of the field to the atom frame. For the scalar (spin 0) field we would have exactly 
the thermal Unruh factor. 

To evaluate the integrals in finite limits, let us suppose again for definiteness that v — lo ^ ^/oloj, = 0, and 
aTe » 1. 

The "counter-rotating" integral Ik{—Lo) does not have stationary points, and its value is 

[v + wY 

The integral Ik{^) is dominated by contribution from the stationary point defined by 

kz . . LO , . 

cosh aTg — smh ar, — — . (48) 

k V 



It is easy to find that 



kz / 'hi ,\A 



k y ayjuj^ - k\c^ 

As in the co-propagating mode case, the ratio (|43|l is anomalously large: it is not exponentially small but linear with 
respect to a/w. 
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VI. COUNTER-RESONANT GAIN AND PARAMETRIC AMPLIFICATION 



Remarkably, not only enhanced spontaneous emission but also laser gain and parametric gain are possible in cavity 
QED via counter-resonant emission by ground state atoms even with random injection times. The gain is reached 
when |/e//ap > 1, or, more exactly, an excess in |/e//(jp over 1 should be greater than the normalized cavity losses. 
For the gain to occur, the time of flight T should be within a certain range to ensure that the atom emits into the 
cavity mode more energy than it takes away, i?2 > Ri- 

In the case of uniformly accelerated atoms, we find that for a co-propagating wave the gain is possible only when 
the acceleration is large enough: a > u. In the opposite case the ratio R2/R1 approaches the asymptotic value 
of a/2TTUj, as was shown in previous sections. Below we plot the ratio R2/R1 for both cases using the incomplete 
gamma- function representation of emission and absorption integrals H29|l . Instead of varying the time of flight T, we 
plot the gain spectrum as a function of the electromagnetic field frequency v for the fixed values of T and the atomic 
frequency lo. 

As is seen from Fig. 2, when oj/a > 1, the emission to absorption ratio drops down to almost zero due to a large 
absorption near the resonance frequency v — oj and then approaches the asymptotic value ajl'KLo in the oscillatory 
way. When ujja < 1, there are strong peaks of a large ratio R2/R1 ^ 1 aX frequencies corresponding to minima of the 
absorption probability; see Fig. 3. Note that the minima of the emission rate are shifted with respect to the minima 
in the absorption. At large field frequencies the envelope of the emission to absorption rate peaks approaches the 
asymptotic value 2{2a/Truj)'^ . 
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FIG. 2; (a) Absorption rate \Ia\^ (thin line) and emission rate |/e |^ (thick line) as functions of the field-to-acceleration frequency 
ratio v/a for the atomic frequency to = 3a and co-propagating wave. Integrals la.e are given by Eq. (11). (b) The ratio of 
emission and absorption rates shown in (a). At large v/a the curve reaches the asymptotic value a/2TTUj. 

The counter-propagating mode is more favorable for the amplification due to sharp dips in the absorption spectrum. 
As is illustrated in Fig. 4, even in the limit uj ^ a the gain spectrum has sharp maxima larger than 1 at the points 
corresponding to nearly vanishing absorption. The case lu < a is qualitatively similar to that of a co-propagating 
mode. 

Note that the peaks of large gain in Figs. 3,4 are not due to maxima of the emission integral but due to minima of 
the absorption probability that are shifted with respect to the minima of the emission spectrum. Absolute values of 
both integrals are small. This is illustrated in the insets to Figs. 3,4 where the emission and absorption spectra are 
shown on the same plot. 

In the optimal regime for amplification, when uj ~ a, v ^ a, and e""-^ <C 1 where the time of flight T ~ L/c, 
one needs to use a longitudinal cavity mode Qn — mrc/L with index n > 1. For example, if aT ~ aL/c — 10, to 
provide v = r2„ = n-Kc/L = 10a one needs n ~ 3. The multimode regime is possible. It is expected to give the same 
qualitative results. 

The effects originated from counter-rotating terms are in fact not uncommon. Two well-known examples are 
parametric resonance and anomalous Doppler effect Q. In all "counter-rotating" processes, an atom can emit a 
photon and simultaneously make a transition from ground to excited state. The required energy is provided by the 
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FIG. 3: (a) Absorption rate \Ia^ (thin line) and emission rate |7e|^ (thick line) as functions of the field-to-acceleration frequency 
ratio v/a for the atomic frequency w = 3a and co-propagating wave. Integrals Ia,e are given by Eq. (11). (b) The ratio of 
emission and absorption rates shown in (a). At large v/a the curve reaches the asymptotic value a/2'Kuj. 
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FIG. 4; (a) Absorption rate \Ia\^ (thin line) and emission rate |/e|^ (thick lino) as functions of the fiold-to-acceleration frequency 
ratio v/a for the atomic frequency ui = 2>a and co-propagating wave. Integrals 7a, e are given by Eq. (11). (b) The ratio of 
emission and absorption rates shown in (a). At large v/a the curve reaches the asymptotic value a/2'Kuj. 



work done by an external force that sustains the center-of-mass motion of an atom along a given trajectory. However, 
an important difference between the nonadiabatic processes considered in this paper and the anomalous Doppler effect 
is that the latter does not require any time-changing parameters. 

It is clear from the above derivation of the omission and absorption probabilities that the enhancement of the 
acceleration radiation is related to a strong nonadiabatic effect at the cavity boundaries. Evidently, this effect 
should exist for an arbitrary trajectory of an atom and in particular, for an atom moving with constant velocity. 
Of course, the presence of acceleration leads to both quahtative and quantitative changes in the excitation rate and 
emission/absorption probabilities by allowing the atom to pass through the resonance between the transition frequency 
of the atom and the Doppler-shifted frequency of the field. 
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For a ground-state atom moving through a cavity with a constant velocity and interacting with a co-propagating 
wave, it is straightforward to obtain the analytic expressions for R2 and 



Ri 



1 



1 



(50) 



i?2 = 



1 



-i{i^'+Lj)T 



= V 



— k • V 
+ k • V 



1/2 



(51) 



Clearly, the factors \/{v' ~ and l/{v' + 10)"^ have the same origin as the nonadiabatic boundary contribution to 
the emission and absorption probabilities given by Eq. (|33l) . When we are far from resonance v' — to, the magnitudes 
of i?i.2 in the cases of constant velocity and constant acceleration are similar and are proportional to the above 
factors. Thus, when the acceleration shifts the frequency ly' further away from the resonance (e.g. when ly < uj ioi the 
co-propagating wave or when u > uj for the counter-propagating wave) , the emission-to-absorption ratio is increasing. 
In this case the effect of acceleration results in the increase of the steady-state number of photons in the cavity as 
compared to the constant velocity case. This tendency is of course reversed when the frequency is shifted towards 
the resonance by acceleration. At the same time, in the case of a constant acceleration we can also have the situation 
when the atom starts far from resonance, then passes through resonance in the course of acceleration, and finally 
ends up far from the resonance. In this case the ratio R2/R1 can be quite large and given by a/2TTUj, while for an 
atom moving with a constant velocity and close to resonance {v' — col <^ ly' , uj the ratio R2/R1 is very small due to a 
strongly enhanced absorption. Thus, depending on the initial conditions, acceleration can lead to either increase or 
decrease in the emission-to-absorption ratio. 

The right-hand side of Eqs. H5U|I . ()51|I strongly depends also on the interference factors e"*'" t^)t g^^g defined 
by the time of flight T, i.e. the phase an atom accumulates relative to the cavity mode while passing through the 
cavity. The ratio R2/R1 can be even greater than one. To achieve R2/R1 > 1, one can tune the time of flight to 
get the proper interference factors: e"'^" -i^)T _^ je^*'-'^ +<^)t _ i| 1, A similar time of flight tuning is used in 
some electronic devices, e.g., in klystrons. The above requirements define a set of the time-of- flight values, with the 
maximum gain corresponding to 



{ly — kv + u})T = (2ni — l)7r; 

{v — kv — Lij)T = 27T,27r, 



(52) 



where ni^2 are integer numbers. For the particular case ni = 0, n2 = — 1 one obtains 2ujT = it. The monochromaticity 
of the beam should satisfy the condition 
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where A — 2t:c/lu, L is a cavity length, and we assumed uj v kv. For v ^ 1 km/s and L ^ X one gets 
Av/v 10^^, which is tough but possible to satisfy. 

The counter-propagating mode is more favorable for the gain since the absorption can then be anomalously small 
while the gain remains as large as for the co-propagating mode. 

Similar interference effects, obviously, are present in the case of a constant acceleration according to Eqs. (|lll) . H29|l . 
(|33|l . and (|34|l . as can be seen in Figs. 1-3. They can also lead to the net gain, as we have already discussed. 

In the case of a parametric resonance, consider an atom moving along an oscillating trajectory z = zq + AcosLUot, 
t = T. The photon absorption and emission probabilities by a ground-state atom ^ are given by 



Jo 

where we introduced a small factor 7 describing the atomic decay. Using 

00 

^ikAcosuot _ Jp[kzA)e^^^°* 

P— — 00 

the above probabilities can be written as 

gJp{k^A) 



(53) 
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(54) 
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where Jp{x) is Bessel's function. Evidently, the probabihties are sharply peaked close to parametric resonance, where 
pLjQ±Lj + Vk — 0- Resonance for emission corresponds to i^+a; = pajg, while the absorption resonance is at v — lo = ptoo- 
When Lo — i>, absorption is always stronger than emission. Indeed, resonance in absorption exists for p — 0, while 
parametric resonance in emission requires p > 1. Therefore, in this case R2/R1 ~ Jpik^A) / jQ{kzA) < 1 for p > 1. 
However, when an atom is not at resonance with the field, one can have parametric resonance in emission but no 
resonance in absorption, which results in the parametric gain. The energy is drawn from the external force causing 
an atom to follow an oscillating trajectory, and the high efficiency of this energy transfer is due to a non-stationary, 
strongly nonadiabatic character of the atomic center-of-mass motion. In the case of Unruh effect, i.e. a uniformly 
accelerated atom in free space, it is also nonadiabaticity that drives simultaneous excitation of the atom and the field. 
However, the efficiency is much lower due to much slower change in the atomic velocity. For an atom entering the 
cavity, a sudden nonadiabatic switch-on of the interaction causes a stronger excitation. 

VII. NONADIABATIC NATURE OF ACCELERATION RADIATION 

The above calculations clearly show that the mechanism of simultaneous excitation of both field and atom is the 
same as for the Unruh effect in free space, namely nonadiabatic transition due to the counter-rotating term ct^ in 
the interaction Hamiltonian lO, i.e. V2. The reason for an enhanced excitation in the cavity is the relatively large 
amplitude for a quantum transition \b,0) — > |a, 1) due to the sudden nonadiabatic switching on of the interaction, 
whereas for the Unruh effect in free space the emission is exponentially small due to a slow switching on. However, 
in both cases there is quite a real emission of a photon accompanied by the excitation of an atom - not just dressing 
of the ground state of an atom as a result of interaction. 

We will now illustrate the above statement by explicit derivation of both the Unruh factor and the enhanced 
excitation factor as a probability of the nonadiabatic transition from the dressed ground state to the dressed excited 
state. Consider first our case of a sudden turn on of the interaction in cavity QED. As a result of the interaction, the 
initial state |6, 0) is no longer an eigenstate of the Hamiltonian. Now, a linear superposition of the excited states of 
the atom and field makes up the dressed ^3 ground state of the interacting system -00 = \b, 0) — -pj^\a, 1) as well as 

the dressed excited state ipi = \a, 1) -I- -p^\b, 0). 

In particular, the amplitude of the bare excited state \a, 1) in ipo is of the order of C ^ iiE' /h{uj + v'). It is easy 
to calculate that the latter corresponds to the atomic excitation probability Paa"" — l^'P ~ \ijlE' /h{uj + ~ l.9^eP, 
where the emission integral le is defined above. This result can be also obtained directly from the density matrix 
equation for the atom, via the atomic counterpart to Eq. (|26|) with a trace over the photon states instead of the 
tratom- This probability has the same origin and value as the well-known Bloch-Siegert shift of a two-level atomic 
transition (13j , Alu/uj — {^.E' /h(uj A- v)Y due to counter-rotating terms in the interaction Hamiltonian. 

The counter-rotating term in Eq. (|33|l represents the contribution from boundaries to the nonadiabatic transition 
amplitudes. In the absence of the boundary contributions, the emission integral Je(w) — Ia{—^) in Eq. H33I) becomes 
exponentially small ^ exp(— Trtj/a) for the small parameter a/2'KUj <C 1 since there are no stationary phase points 
in the integration interval. The absorption integral la does have a point of stationary phase when the atomic 
frequency uj is brought into resonance with the field due to the time-dependent Doppler shift of the mode frequency 
[1^ v'{t) = i^exp(— ar). This fact explains why the related exponential factor effectively disappears from the 
absorption integral H34() when a <C 2-ku). As a result, if there are no edge effects, we obtain the same excitation factor 
R2/R1 ~ exp(— 27rijj/Q!) as in the Unruh effect (in free space). This means that in order to observe the standard 
Unruh result one has to extend the mode profile g{z) near the boundaries, i.e., eliminate nonadiabatic boundary 
contributions. 

Similary to what we did for the sudden turn-on case, let us now demonstrate the nonadiabatic nature of the Unruh 
effect by the following explicit derivation of the Unruh factor as a probability of the nonadiabatic transition t/jq V'l 
from the dressed ground state. The Shroedinger equation ihdij^/dT = Hip in the two-level case ip = cqV'o + ciipi 
yields dci/dr + {iEi/h+ {'tpi\'tpi))ci = — co(^olV'i)- The difference between the eigenenergies is, to the first order. 
El — Eq ~ h{Lu + ly'). For small nonadiabatic coupling — (■0o|'0i) = ( J+^J' ) ^ to + v' , the perturbation solution 

is |cip — I /J" exp[i {v' + u))dT"]-^ (5+1^) '^''"'P- If n^'^ make the assumption of an adiabatic switching (on 
and off) of the interaction (/(r) as in standard Unruh effect treatments, then after integration by parts the latter 
integral is reduced to the integral Ie{uj) = Ia{—oj) in Eqs. (6) but in the infinite limits, i.e. without edge effects. This 
yields the standard Unruh factor |cip cx exp(— 27ra;/a). This derivation clearly shows the dramatic effect of boundary 
contributions leading to a large amplitude ^ g(T)/{uj + v') of the atomic excited state |a). Only if we eliminate the 
edge effects by adiabatic switching of the interaction do we retrieve the exponentially small excitation factor. 

Note that in the cavity the excitation factor exp(— fizz/fesT) = R2/R1 — a/2'KUj is determined by the first power 
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of the same nonadiabaticity parameter a/2TTij. The reason for this effect is the existence of a true resonance, i.e., a 
stationary-phase point, in the absorption coefficient. As mentioned earlier, this yields a resonance between the atomic 
transition frequency and the Doppler-shifted frequency of the field seen by the atom, u) + ^(^e"""^) ~ 0, and is 
responsible for the aforementioned effect. 



VIII. CONCLUSIONS 

Our simple model clearly demonstrates that the ground state atoms accelerated through a vacuum-state cavity 
radiate real photons. For relatively small acceleration a < 2nujc. the excitation Boltzman factor exp(— /tiz/fesT) ~ 
a/2'KU! is much larger than the standard Unruh factor exp(— 27rci-'/Q). The physical origin of the field energy in the 
cavity and of the real internal energy in the atom is, of course, the work done by an external force driving the centcr-of- 
mass motion of the atom against the radiation reaction force. Both the present effect (in a cavity) and standard Unruh 
effect (in free space) originate from the transition of the ground state atom to the excited state with simultaneous 
emission of photon due to the counter-rotating term a~^^~^ in the time-dependent Hamiltonian (9). Thus, these effects 
have essentially the same counter- resonant, nonadiabatic mechanism. We emphasize that there is emission of real 
photons in both cases; however the emission probability is exponentially small for the standard Unruh condition of the 
absence of boundaries and slow turn-on of the interaction - simply because the nonadiabatic effect is very small in the 
latter case. The enhanced rate of emission into the cavity mode comes from the enhanced nonadiabatic transition at 
the cavity boundaries; the standard Unruh excitation comes from the nonadiabatic transition in free space due to the 
time dependence of the Doppler-shifted field frequency v' = z^e"""^, as seen by the atom in the course of acceleration. 
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Figure Captions 

Fig. 1. (a) Atoms or ions in their ground state \b) are accelerated through a single- mode microwave or optical 
cavity, (b) Resonant absorption or emission: an atom is excited (deexcited) as it simultaneously absorbs (emits) 
a photon, (c) Counter-resonant absorption or emission processes that are usually neglected in the "rotating wave 
approximation" : an atom is excited (deexcited) as it simultaneously emits (absorbs) a photon, (d) the energy for 
counter-resonant emission is drawn from work done by a force accelerating an atom. 

Fig. 2. (a) Absorption rate |/ap (thin line) and emission rate |/ep (thick line) as functions of the field-to- 
acceleration frequency ratio v/a for the atomic frequency w = 3a and co-propagating wave. Integrals Ia,e are given 
by Eq. (11). (b) The ratio of emission and absorption rates shown in (a). At large v/a the curve reaches the 
asymptotic value a/ 2^0). 

Fig. 3. (a) Absorption rate (thin line) and emission rate |/ep (thick line) as functions of the field-to- 

acceleration frequency ratio f/a for the atomic frequency lv = a/3 and co-propagating wave. Integrals Ia,e are given 
by Eq. (11). The inset shows the tails in more detail, (b) The ratio of emission and absorption rates shown in (a). 
At large v/a the curve reaches the asymptotic value (2a/7ra;)^/2 = 8. 

Fig. 4. (a) Absorption rate |7aP (thin line) and emission rate |7eP (thick line) as functions of the field- to- 

acccleration frequency ratio v/a for the atomic frequency lo = 3a and counter-propagating wave. Integrals /(,,(. arc 
given by Eq. (11). The inset shows the tails in more detail, (b) The ratio of emission and absorption rates shown in (a). 



